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Zinc  oxide  varistors^  are  made  of  ceramic  materials  consisting  of  zinc  oxide 
with  a  few  percent  of  additive  oxides,  mostly  bismuth  or  praseodymium  oxide 
(varistor-forming  ingredient)  and  transition  metal  oxides 
(performance-enhancing  ingredients).  The  varistor  shows  a  strongly  nonlinear 
i-V  characteristics,  with  a  sudden  transition  from  insulating  to  conducting 
state  at  the  "breakdown"  voltage.  This  sharp  transitioh' makes  it  an  ideal 
device  to  protect  electrical  equipment  against  sudden  voltage  surge  or 
electromagnetic  pulse. 


(mS) 


Figure  1 

When  a  square  voltage  pulse  of  the  mplitude  greater  than  the  breakdown 
voltage  is  applied,  the  varistor  curre.  t  shows  a  capacitance  charging  peak 
followed  by  a  minimum  before  it  rises  to  the  steady  state  value  as  in  Fig.  1 . 

The  current  rise  time  is  usually  of  the  order  of  1  ps  except  at  the  applied 
voltage  very  close  to  the  breakdown  value;  above  a  certain  voltage  for  a  given 
sample,  the  rise  time  remains  constant  and  independent  of  the  applied  bias 
This  current  minimum  or  voltage  overshoot  phenomenon  is  found  for  all  ZnO 
varistors  with  highly  nonlinear  i-V  characteristics.  A  few  parts  per  million 
of  indium  removes  the  voltage  overshoot  but  it  also  causes  a  drastic  reduction 
of  the  nonlinearity^. 

The  voltage  overshoot  effectively  increases  the  turn-on  time.  In  order  to 
design  a  varistor  with  a  fast  turn-on/turn-off  time  and  high  nonlinearity,  a 
good  understanding  of  microscopic  mechanisms  for  the  voltage  overshoot  and 
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high  nonlinearity  becomes  essential.  This  grant  supported  a  preliminary 
investigation  of  the  varistor  transient  behaviors  associated  with  turn-on  and 
turn-off  of  the  varistor. 


of  the  most  important  results 


We  have  studied  the  temporal  response  of  varistor  devices  fabricated  by  GE  as 
well  as  Toshiba  varistors.  We  have  comfirmed  the  tim'e  responses  reported  by 
F.  Modine  of  Oak  Ridge  and  further  explored  the  initial  fast  decay  region,  Fig.  2. 
We  find  that  after  the  initial  exponential  decay  there  is  a  power  law  behavior  / 
where  n  is  close  to  unity  for  the  applied  bias  below  the  breakdown 
(threshold)  voltage,  very  similar  to  the  transient  current  following  the  removal 
of  the  bias  voltage^. 


Figure  2 


The  temperature  dependence  of  the  transient  current  is  small-  too  small  to  be 
associated  with  charge  carrier  creation  mechanisms.  In  stead  of  one 
activation  energy  reported  in  the  past,  a  new  analysis  of  the  Modine's  data 
yields  two  activation  energies  of  ~  10  meV  and  ~  160  meV  for  the  delayed 
transient  current  in  the  region  from  10  second  to  10^  second.  In  the  region 
immediately  following  the  exponential  decay,  our  own  data  show  an  activation 
energy  of  ~  160  meV;  the  lower  activation  energy,  ~  10  moV,  has  not  been  seen 
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because  we  have  not  been  able  to  carry  out  measurements  at  low  temperatures 
where  it  is  dominant  over  the  higher  activation  energy.  These  activation 
energies  are  too  small  to  be  consistant  with  the  slow  decaying  nature  of  this 
polarization  current  if  they  are  associated  with  carrier  creation  processes; 
they  must  be  associated  with  the  carrier  mobility.  We  have  developed  a  theory 
which  explains  the  slow  decaying  polarization  current.  We  propose  that  this 
slow  decaying  polarization  current  is  generated  by  the  tunnelling  into  the 
conduction  band  of  deep  donors  under  the  influence  of  the  depletion  layer 
potential  energy  provided  by  the  shallow  donors  which 'ionize  in  a  very  short 
time,  say  a  few  nanoseconds,  when  the  Fermi  energy  is  shifted  down  across 
them.  Since  this  is  only  a  part  of  our  theoretical  program,  we  start  with  a 
general  description  of  our  view  of  the  varistor  mechanism. 

Since  existing  theories  do  not  provide  satisfactory  explanations  for  the  long 
turn-on  (the  voltage  overshoot)  and  short  turn-off  time,  and  also  how  some 
transition  metal  additives  (performance  ingredients)  improve  the  varistor 
performance  so  drastically  by  increasing  the  nonlinearity,  we  have  tried  to 
develop  a  theory  which  provides  explanations  for  them.  The  theory  we  are 
developing  is  based  on  the  following  observation: 

1 .  The  double  Schottky  barrier  model  of  the  grain  boundary  region 
is  valid, 

2.  The  performance  ingredients  such  as  Mn  and  Co  introduce 
additional  electron  donor  levels  in  the  energy  band  gap, 

3.  There  are  shallow  donors  within  a  few  kT  of  the  conduction  band 
edge, 

4.  Tunnelling  as  well  as  thermally  activated  electronic  processes 
are  involved,  and 

5.  Hot  electrons  have  short  mean  free  path  and  consequently  a  short 
thermalization  length. 

We  start  with  the  double  Schottky  barrier  models  of  Pike^  ;  Mahan,  Levinson 
and  Philipp®:  and  Blatter  and  Greuter^,  but  recognize  that  the  short  mean  free 
path,  with  accompanying  energy  loss,  make  the  creation  of  holes  in  the  valence 
band  --  or  impact  ionization  -  by  the  electrons  in  the  conduction  band  unlikely 
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at  the  bias  of  approximately  3  eV  if  the  width  of  the  depletion  layer  remains  as 
large  as  1000  A  or  greater.  We  estimate  the  mean  free  path  as  small  as  10  A, 
leading  to  the  conclusion  that  at  least  100  scattering  events  occur  before  the 
electron  leaves  the  depletion  layer.  The  electrons  coming  over  the  barrier  may 
be  assumed  thermalized  by  the  time  they  leave  the  depletion  layer.  We 
consider  the  possibility  that  the  width  of  the  depletion  layer  car  decrease  on 
the  positively  biased  side  as  the  deep  donors  ionize  by  tunneling  into  the 
conduction  band  and  also  by  impact  ionization  by  hot  electrons.  Although  there 
are  donors  at  various  energy  levels,  we  assume  only  two  types  of  donors, 
shallow  donors  of  density  n-|  and  deep  ones  of  n2,  in  order  to  bring  out  the 

physics  clearer. 


Figure  3 


With  zero  bias,  all  donors  above  the  Fermi  level  are  ionized  as  illustrated  in 
Figure  3.  When  a  bias  V  is  applied,  the  Fermi  level  on  the  positive  biased  side 
is  lowered  while  that  of  the  negative  side  is  raised.  Soon  after  the  bias  is 
applied,  some  conduction  electrons  move  into  the  depletion  layer  on  the 
negative  side  and  neutralize  some  ionized  donors  while  the  shallow  donors  on 
the  positive  side  which  are  now  above  the  Fermi  level  become  ionized  in  a  time 
of  x-j  while  the  deep  donors  will  take  a  longer  time  X2  to  ionize.  Soon  after, 

that  is,  x-|  after  the  bias  is  applied,  the  conduction  band  edge  in  the  depletion 
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layer  (the  electron  potential  energy  curve)  consists  of  two  parts;  the  sum  of 
the  parabolic  curve  due  to  n-|  and  the  other  parabola  due  to  r\2  as  given  by  Eq.  2. 

(e2ni/2e)  (X  -  xpj)2  fo^xp^  <  x  <xp 

V(x)  -  Vp  =  (2) 

(e2ni/2e)  (x  -  xp)^  +  (e2n2/2e)  (x  -  xp2o  for  0  <  x  <  xp2o 

where  Vp  is  the  barrier  height  and  xp2o  is  the  width  at  zero  bias  of  the 
ionized  deep  donor  distribution,  for  the  positive  side. 

All  the  shallow  donors  above  the  positive  side  Fermi  level  are  ionized  to 
provide  a  wide  parabola  while  only  ionized  deep  donors  are  those  which 
were  ionized  at  zero  bias  and  provide  a  parabola  of  the  zero  bias  width,  Eq. 

3. 

XRo  =  [  2eV0/e2{n^  +  02)]"'^^  (3) 

where  Vp  is  the  barrier  height  at  zero  bias. 

After  X  2  the  deep  donors  above  the  positive  side  Fermi  level  would  have 

been  ionized,  resulting  in  the  narrowing  of  the  depletion  layer  as  indicated 
in  Figure  4.  If  n2  were  much  larger  than  n-i ,  the  narrowing  would  be 

drastic  and  tunneling  from  the  grain  boundary  may  become  feasible. 
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Figure  4a 


Figure  4b 


As  the  bias  is  increased,  some  electrons  coming  over  the  barrier  and 
across  the  grain  boundary  may  have  enough  energy  to  excite  a  valence  band 
electron  to  the  conduction  band  creating  a  hole.  This  could  happen  for  the 
bias  comparable  to  the  band  gap  only  if  the  depletion  layer  width  is  small 
enough  so  that  the  electron  does  not  lose  significant  fraction  of  its  energy 
before  it  leaves  the  depletion  layer.  The  value  for  n-)  is  often  equated  to 

the  estimated  carrier  density  of  1 0^ ^  cm'^  .  If  the  value  of  n2  is  1 0^® 

cm'3,  the  depletion  layer  width  would  shrink,  in  X2,  to  I/V1OOO  «  1/32 

of  the  width  at  x-j .  A  1 000  A  wide  depletion  layer  would  shrink  to  a  31  A 

wide  thin  layer,  making  it  feasible  for  an  electron  move  across  the  layer 
without  losing  more  than  a  fraction  of  its  energy  and  thus  making  hole 
creation  possible  when  the  barrier  height  on  the  positive  side  is  slightly 
larger  than  the  band  gap.  The  value  of  10^®  cm'3  for  n2  is  not 

unreasonably  large  since  a  few  percent  of  transition  metal  oxides  are 
added  to  improve  the  varistor  performance.  This  mechanism  can  provide 
explanations  to  (i)  why  the  turn-on  voltage  per  grain  boundary  (the 
breakdown  voltage)  seems  to  be  independent  of  the  additives,  (ii)  why  the 
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metai  oxide  varistor  with  high  nonlinearity  coefficient  has  a  voltage 
overshoot  lasting  microseconds,  and  (iii)  why  the  displacement  current 
peak  is  followed  by  a  very  small  transient  current  which  decays  over  a 
long  time  if  the  bias  is  less  than  the  turn-on  voltage  and  also  when  the 
applied  bias  is  turned  off. 

As  stated  previously  a  slow-decaying  transient  current  follow  the  removal 
of  the  bias  and  also  follow  the  displacement  current  peak  for  the  bias  less 
than  the  breakdown  (turn-on)  voltage.  This  transient  current  has  been 
found,  by  Modine  and  us,  to  have  a  power  law  time  dependence  /  «:  1/f 
from  a  few  microseconds  to  thousands  of  seconds,  with  the  value  of  n 
close  to  unity.  The  power  law  time  dependence  can  be  explained  if  we 
assume  that  the  charge  carriers  are  generated  by  the  tunnelling  into  the 
conduction  band  of  the  deep  donor  electrons  which  are  now  metastable 
since  the  Fermi  level  shifted  down  across  their  energy  levels.  As 
discussed  above,  the  deep  donors,  after  _  are  under  the  influence  of  the 

parabolic  potential  provided  by  the  ionized  shallow  donors.  We 
approximate  the  potential  barrier  by  a  triangle  (constant  potential 
gradient  approximation)  as  shown  in  Fig.  7  and  Eq.  4. 

V(^)  -  E  =  [Ed/(X2  -  X)]  (X2  -  k)  (4) 

for  the  trap  at  x.  E<j  =  V(x)  -  E  is  the  trap  depth  measured  from  the 
conduction  band  edge  at  x.  The  WKB  approximation  yields 

X  oc  exp  { [(2mE(j/h2]''/2(.4££^/3g2p,^x)}  (5) 

where  x  takes  negative  values  because  of  our  choice  of  the  zero.  An 
interesting  consequence  of  this  formula  is;  the  life-time  of  the  trap  is 
greater  for  those  farther  inside  the  bulk,  i.e.  for  smaller  (-x).  This 
dependence  of  the  life-time  on  x  does  produce  the  power  law  time 
dependence  of  the  slow-decaying  transient  current  in  a  time  range. 

We  may  assume  the  current  to  be  proportional  to  the  rate  of  carrier 
creation  or  carrier  injection  from  the  depletion  layer  to  the  bulk.  Then 
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where  is  the  life-time  at  x  =  adjacent  to  the  grain  boundary. 
From  Eq.  5,  we  derive 


dx  _ 
dt  ^ 


1  1 


2 


_1_ 


X 


oc 


2^ 

-r 


(7) 


where  A  is  a  constant  independent  of  x ;  the  weak  t-dependence  in  the 
brackets  is  neglected.  Substituting  Eq.  7  into  Eq.  6  yields 


i  it) 


—  dx  =  —  ( 1  -  e 

X  t 


-t/t 


for  t  »  Xb 


(8) 


A  slightly  different  formulation  is  given  in  our  paper  ^  to  be  published  in 
Advances  in  ceramics. 
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